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ABSTRACT 

Over the last few years several models have been proposed to interpret the widespread 
soft excess observed in the X-ray spectra of type 1 active galactic nuclei (AGN). In 
particular, reflection from the photoionized accretion disc blurred by relativistic effects 
has proven to be successful in reproducing both the spectral shape and the variability 
pattern of many sources. As a further test to this scenario we present the analysis 
of a recent ^100 ks long Suzaku observation of Arakelian 120, a prototypical 'bare' 
Seyfert 1 galaxy in which no complex absorption system is expected to mimic a soft 
excess or mask the intrinsic properties of this key component. We show that a re- 
flection model allowing for both warm/blurred and cold/distant reprocessing provides 
a self-consistent and convincing interpretation of the broadband X-ray emission of 
Ark 120, also characterized by a structured iron feature and a high-energy hump. 
Although warm absorbers, winds/outflows and multiple Comptonizing regions may 
play significant roles in sources with more spectral complexity, this case study adds 
evidence to the presence of blurred disc reflection as a basic component of the X-ray 
spectra of type 1 AGN. 
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1 INTRODUCTION 

Over the 2-10 keV energy range, the X-ray spectra of type 1 
active galactic nuclei (AGN) are usually well-reproduced by 
a power-law component, which is believed to originate from 
the Comptonization of cold seed photons in a hot coronal re- 
gion above the accretion disc, according to the so-called two- 
phase model (Haardt & Maraschi 1991, 1993). When extrap- 
olated to the lower energies, the simple power-law trend fails 
to account for the extra emission systematically observed in 
unabsorbed Seyfert galaxies (after Arnaud et al. 1985). Due 
to its extremely smooth shape, this soft excess turns out to 
be consistent with several physical models. In analogy with 
the high/soft states of galactic black hole binaries (Remil- 
lard & McClintock 2006), dominated by thermal emission 
from the disc, the initial attempts at modelling the AGN 
soft excess were based on single or multicolour blackbody 
components. Although well-fitting, this thermal continuum 
requires a temperature far higher than that expected for a 
standard thin disc around a supermassive black hole. More- 
over, when tested on large AGN samples, the characteristic 
temperature is found to be remarkably constant over almost 
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four orders of magnitude in the black hole mass, in spite of 
very different Eddington rates (e.g. Crummy et al. 2006, and 
references therein; Miniutti et al. 2009). A thermal soft ex- 
cess does not even comply with the Stefan-Boltzmann law 
for the largest variations in the overall luminosity (Ponti 
et al. 2006), and the inferred size of the emitting region is 
usually unreasonably small with respect to the gravitational 
radius (r g = GM/c 2 ). 

In terms of spectral decomposition, similar results can be 
obtained by invoking the Comptonization of disc photons 
in a cold and optically thick plasma (Page et al. 2004). 
This raises questions on the nature of the X-ray corona. 
In the case of two detached Comptonizing regions with dif- 
ferent physical properties (e.g. Dewangan et al. 2007), most 
of the limitations of the thermal scenario described above 
also apply to the 'soft' corona. On the other hand, a single 
plasma with hybrid thermal/non-thermal electron distribu- 
tion (Coppi 1999) may give rise to both the soft excess and 
the hard power-law component. 

Alternatively, the universal shape of the soft X-ray emission 
has been linked to atomic physics: Gierlinski & Done (2004) 
suggested that the soft excess is actually a fake continuum 
component due to a broad absorption trough at ~2-5 keV, 
arising in partially ionized gas along the line of sight subject 
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to high velocity smearing. The latest simulations, however, 
prove that the properties of any realistic accretion disc wind 
are not able to reproduce the observed smoothness of the 
soft excess (Schurch & Done 2007; Schurch, Done & Proga 
2009). Another viable explanation is that of reflection from 
the photoionized surface layers of the accretion disc itself, 
where the relativistic motion of the infalling matter pro- 
vides the blurring of the narrow atomic features (Fabian 
et al. 2002; Crummy et al. 2006). This model occasionally 
implies a strong suppression of the intrinsic power law in 
order to account for a prominent soft excess, as justified in 
the context of strong gravitational light bending (Miniutti 
& Fabian 2004). 

On sheer statistical grounds all the interpretations outlined 
so far yield acceptable results (Sobolewska & Done 2007). 
To pursue further this kind of study, access is needed to 
energies beyond ~10 keV, where the physical models make 
different predictions. The advent of Suzaku (Mitsuda et al. 
2007) and its hard X-ray detector (HXD; Takahashi et al. 
2007) has made it possible to put solid spectral constraints 
up to ~70 keV, opening a new era of AGN observations. 
Concerning the soft excess, it should also be noticed that 
the possible presence of complex and variable absorption can 
mask the intrinsic appearance of this critical component. In- 
deed, warm absorbers are very common among type 1 AGN 
(Crenshaw, Kraemer & George 2003; Blustin et al. 2005). 
It is therefore desirable to select a target with the cleanest 
view of the nuclear regions. 

Arakelian 120 (z = 0.0327) is a rare case of a 'bare' Seyfert 
galaxy, hence it represents the optimal candidate to test 
whether the soft excess is a signature of blurred reflection 
from the inner accretion disc. No evidence for reddening is 
found in the infrared (Ward et al. 1987), and ultraviolet 
observations establish that Ark 120 is devoid of intrinsic ab- 
sorption (Crenshaw et al. 1999). In the soft X-ray band the 
source has been shown to have a steep spectrum by both 
EXOSAT (Turner & Pounds 1989) and ROSAT (Brandt et 
al. 1993), while the XMM-Newton Reflection Grating Spec- 
trometer data allow stringent upper limits (~l-2 orders of 
magnitude lower than those of usual Seyfert Is) to be placed 
on the ionic column densities of any possible warm absorber 
(Vaughan et al. 2004). It is worth emphasizing that Ark 120 
is a broad-line Seyfert 1 (BLS1) galaxy, whereas the soft ex- 
cess has been long associated with narrow-line sources only. 
NLSls are rather eccentric objects, known for their very 
steep X-ray spectrum (Boiler, Brandt & Fink 1996), high 
accretion rate (Grupe 2004, and references therein), large- 
amplitude X-ray variability on short timescales (Gallo et 
al. 2004), metal overabundance (Shemmer & Netzer 2002; 
Fabian et al. 2009) and enhanced star formation (Sani et 
al. 2010). Ark 120 is an outstanding counter example of a 
normal BLS1 with a prominent soft excess, ensuring as such 
a completely unbiased exploration of this component. 
This work is organized as follows: Section 2 concerns the ob- 
servation and data reduction; our results are presented and 
fully discussed in Section 3, while in Section 4 we summarize 
our conclusion and outline the future research. 



2 OBSERVATION AND DATA REDUCTION 

Ark 120 was observed by Suzaku on 2007 April 1-3 in the 
HXD nominal position, with a resulting net exposure of 
~101 ks for the X-ray imaging spectrometer (XIS; Koyama 
et al. 2007) and ~89 ks for the HXD/PIN detector. Events 
were collected by the three operational XIS CCDs in both 
3x3 and 5x5 editing modes (71 and 30 ks exposures, re- 
spectively). Following the standard procedure illustrated in 
the Suzaku Data Reduction Guidejj we used the HEASOFT 
FTOOLS 6.8 package to obtain a new set of clean event files 
for each detector, editing mode and telemetry by reprocess- 
ing the unfiltered events with the latest calibrations. The 
source spectra and light curves were extracted from circular 
regions with radius of ~3.5 arcmin (200 pixels) centred on 
the target, while the background was evaluated on the same 
chip from adjacent regions devoid of significant contamina- 
tion. Finally, the source and background spectra from the 
two front-illuminated detectors (XIS0 and XIS3), as well as 
the response files generated through the 'xisresp' script, were 
merged and rebinned by a factor of 4. 

Concerning the HXD/PIN data reduction, we again repro- 
cessed the unfiltered event files using the standard tools and 
got the output spectrum by running the 'hxdpinxbpi' script, 
which takes into account the contribution of both the non 
X-ray and the cosmic X-ray background and applies the 
dead time correction. This returned a source count rate of 
(12.1 ± 0.3) x 10~ 2 s" 1 , corresponding to 18.7 per cent of 
the PIN total counts. 

The spectral analysis has been performed using the xspec 
vl2.6 fitting package, and involves only the 0.5-12 keV en- 
ergy range of the two front-illuminated XIS detectors; in- 
deed, the 1.7-2.0 keV interval appears to be strongly af- 
fected by systematic calibration uncertainties around the 
instrumental silicon K-edge, hence it has been excluded. 
The XIS1 back- illuminated spectrum has been employed 
throughout as an independent check. The source is confi- 
dently detected at high energy up to ~40-50 keV; we have 
therefore considered the conservative 12-40 keV range of 
the HXD/PIN spectrum. In order to allow the use of \ 2 
minimization during the spectral fitting, the XIS data were 
grouped so that each energy channel contains no less than 20 
counts, while a minimum of 50 counts per bin was adopted 
for the HXD/PIN spectrum. The uncertainties reported in 
this work correspond to the 90 per cent confidence intervals 
(A^ 2 = 2.71) for the single parameter of interest. Fluxes 
and count rates are given at the la level. 



3 DATA ANALYSIS AND DISCUSSION 

In Fig.[T]we show the data/model ratio plot obtained by fit- 
ting the Suzaku spectrum of Ark 120 with a simple power law 
over the 2.5-5.5 keV range. The output is obviously model- 
dependent, especially with respect to the magnitude of the 
soft and hard excess; none the less, both the sizable spec- 
tral curvature and the prominent iron emission are clearly 
brought out. Our analysis is first aimed at a detailed inspec- 
tion of the latter feature, and is subsequently extended to 
the whole 0.5-40 keV energy range. Also, we briefly examine 
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Figure 1. The Suzaku XIS0+XIS3 and PIN spectra of Ark 120 
plotted as a data/model ratio after fitting the 2.5-5.5 keV range 
only with a power law. This sketch, though model-dependent, is 
straightforward in bringing out both the prominent iron feature 
at ~6 keV and the overall spectral curvature, the latter suggest- 
ing the presence of a smooth soft excess below 2 keV and of a 
broad reflection hump beyond 10 keV. (The spectral data shown 
throughout this work have been rebinned for plotting purposes 
only). 

the source variability, and review the previous XMM-Newton 
observation of Ark 120 in the light of the new evidence from 
the present work. 

3.1 Iron K-shell emission 

After excluding the 5-7.5 keV region characterized by the 
complex iron feature, the 2.5-12 keV spectrum has been fit- 
ted with a basic model consisting of a power law modified 
by Galactic absorption only. Despite the apparent lack of 
local obscuration in Ark 120, the foreground column den- 
sity towards the source is quite large due to its low Galactic 
latitude (b = — 21?13). Two different values for A^h have 
been determined from the A21-cm neutral hydrogen emis- 
sion maps, i.e. 0.98 x 10~ 21 cm -2 (Kalberla et al. 2005) and 
1.26 x 10~ 21 cm" 2 (Dickey & Lockman 1990). Throughout 
this work, Galactic absorption has been modelled by means 
of the TBABS code (Wilms, Allen & McCray, 2000), assum- 
ing the estimates above to be the extremes of the allowed 
range for Ah- We will briefly return to this point later on, 
since for now the selection of a precise value of A^h does not 
affect the study of iron emission around 6 keV. 
The shape of this feature is shown in the top panel of Fig. [5] 
Different attempts have been made in order to disentangle 
the components of such a structured profile. Here we con- 
sider the following four cases: two narrow lines; three narrow 
lines; two narrow lines and a disc line; two broad lines. 
1) A model including only a couple of narrow gaussian lines 
of width a = 10 eV yields a statistically good fit to the 
data, with x^/d.o.f. = 0.869/283. The lines are centred at 
~6.42 and 6.96 keV in the source frame, and are consistent 
with the energy expected for K-shell emission from neutral 
and H-like iron; their equivalent widths (EW) are ~90 and 
40 eV, respectively. However, this solution is clearly insuffi- 
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Figure 2. Top panel: the structured profile of iron K-shell emis- 
sion obtained as a data/model ratio after fitting the 2.5-5 plus 
7-12 keV energy range with a single power law. It is evident that 
two narrow gaussian components are not sufficient to reproduce 
the whole feature. The data from XIS0 (blue dots), XIS1 (red 
squares) and XIS3 (green crosses) are shown separately. Bottom 
panel: same as above, after removing the disc line component 
from the best fitting-model. 



cient to reproduce the observed feature. 

2) By adding another narrow line, the fit undergoes a consid- 
erable improvement (x 2 /d.o.f. = 0.789/281). The new line 
has an energy of ~6.67 keV, driving a straightforward iden- 
tification as He-like iron emission, and an EW of ~30 eV. 
But also this interpretation is not fully adequate. 

3) In fact, the most successful way to fit the residuals left by 
the two narrow lines is to allow for a disc line, whose broad- 
ening is due to relativistic effects. By including in the model 
a LAOR profile (Laor 1991) we get an excellent fit of the whole 
iron emission feature (x 2 /d.o.f. = 0.686/280). The shape 
of this skewed component is shown in the bottom panel of 
Fig. and clearly exhibits the red elongated wing typical 
of lines arising from the inner disc, even if not as prominent 
as in the well-known case of MCG-6-30-15 (Miniutti et al. 
2007). For this reason, it is not possible to constrain simul- 
taneously all the line parameters (see Table[T|: these include 
the disc inner radius (n n ) and inclination with respect to the 
line of sight and the emissivity index q (under the assump- 
tion of a radial emissivity profile e(r) cx r~ q ). None the less, 
we note that the best-fitting inner radius of ~13 r g does 
not necessarily require the extreme gravity regime typical of 
rapidly rotating black holes. It is important to stress that 
alternative explanations invoking a spectral upturn due to 
complex absorption effects, such as those proposed for the 
same MCG-6-30-15, are definitely not viable in the case of 
Ark 120. 

4) In the latter stage, the properties of the two narrow 
lines are exactly the same as found in the original model. 
Thus, it is worth investigating another possibility by re- 
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Table 1. Best-fitting parameters for different models of the iron K-shell feature: 1) two narrow lines; 2) three narrow lines; 3) two narrow 
lines plus a LAOR disc line; 4) two broad lines. The asterisk-marked parameters have not been allowed to vary. 
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taining the basic template and allowing the width of the 
two lines to vary. The consequent fit refinement is not 
as significant as when a disc line is involved, being now 
xl/d.o.i. = 0.735/282: the difference of A\ 2 ^ -15.2 ob- 
tained with the loss of two degrees of freedom is not likely to 
be simply a chance improvement and supports the relativis- 
tic line detection. Moreover, although the double-peaked en- 
ergy of the blended feature is still consistent within the er- 
rors with neutral and H-like iron emission (Table [TJ, the 
resulting width of a — 113(±21) eV poses the question 
about the physical location wherein these lines arise. Such a 
value, in fact, corresponds to a full width at half-maximum 
(FWHM) broadening of ~ 12 x 10 3 km s" 1 , which is a factor 
of ~2 larger than that observed in the optical permitted lines 
(FWHM H/3 ~ 5800 km s" 1 ; Wandel, Peterson & Malkan 
1999). This could hint at a sort of X-ray broad-line region 
(BLR) internal to the optical one. 

Any further discussion on the possible origin of these lines 
(either broad or narrow) is deferred to the next section, in 
which we address this issue within the context of X-ray re- 
flection models. We are confident that relativistic effects are 
in place and have to be taken into account also when consid- 
ering the entire spectral range, since all the interpretations 
of the iron emission profile in Ark 120 making no resort to 
a disc line turn out to be less successful and convincing. 

3.2 Blurred reflection model 

We now extend our analysis to the whole 0.5-40 keV en- 
ergy range, in order to understand the origin of the spectral 
curvatureQ Excess emission beyond ~20-30 keV is usually 

2 The cross calibration between the XIS and PIN spectra has 
been fixed to the recommended value of 1.17. 



interpreted as due to the reprocessing of the primary X-ray 
radiation: the combination of photoelectric absorption and 
Compton scattering in the illuminated material gives rise to 
a broad reflection hump (e.g. George & Fabian 1991, and 
references therein). Besides this additional continuum com- 
ponent and iron fluorescence, below ~2 keV the reflected 
spectrum is expected to be dominated by a wealth of emis- 
sion lines from oxygen and other abundant elements, like C, 
N, Ne, Mg, Si, S (Ross & Fabian 1993). Since it is fairly con- 
ceivable that in many cases the accretion flow itself acts as 
the most efficient mirror, depending on the ionization stage 
of the disc outer layers and on the relativistic motions of the 
inner regions, the stack of individual features can be blurred 
into the smooth shape of the soft excess. 
As stated above, the slight spectral curvature observed at 
~6 keV in Ark 120 cannot be explained as the product 
of (multiple) covering effects, suggesting instead the pres- 
ence of a broad skewed profile within the iron emission fea- 
ture and indicating a strong gravity regime. It is therefore 
reasonable to include in our general model two reflection 
components, for which we have used the self-consistent re- 
FLIONX table models of Ross & Fabian (2005): the first one 
is expected to arise from almost neutral material at great 
distance from the X-ray emitting region, likely situated on 
the dusty torus scale. The second one can be ascribed to the 
partially ionized surface of the disc, and has been convolved 
with the kdblur kernel in xspec to account for relativistic 
effects. The model obviously comprises the primary power 
law, absorbed by the Galactic column density, and also the 
two unresolved gaussian lines identified above, correspond- 
ing to neutral and H-like iron emission. Galactic absorption 
is found to be consistent with Ah = 0.98 x 10 -21 cm -2 , and 
has been frozen accordingly. 

The reflection templates include among their free parame- 
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Figure 3. Top panel: Suzaku XIS and HXD/PIN spectrum of 
Ark 120 and best-fitting reflection model (light green solid line). 
The contribution of all the different components is disentangled: 
primary power law (blue), blurred (red, dashed) and distant (deep 
green, dot-dashed) reflection, narrow iron lines (magenta). The 
shaded 1.7-2 keV region has been excluded from the fit. Bottom 
panel: residuals in units of a (the error bars have size one). 

ters iron abundance (which has been assumed to be single- 
valued across the source) , ionization state (described by the 
quantity £ = 4n F/rm, where Tin is the hydrogen number 
density of the gas exposed to an X-ray flux F) and photon 
index of the illuminating power law (which has been taken 
to be one and the same with the direct power law). The 
key blurring parameters are those already introduced above 
for the LAOR shape, i.e. inner radius n n , emissivity index q 
and disc inclination. The disc outer radius is usually hard 
to constrain and has been fixed to 400 r g . The normaliza- 
tions of the power-law and reflection components are not 
tied to each other: in fact, in the presence of light bending 
effects (e.g. Miniutti & Fabian 2004) which may focus the 
illuminating radiation on to the disc violating isotropy, the 
traditional reflection fraction has no clear link with the geo- 
metrical covering factor. As a consequence, in this work we 
simply weight the contribution of the reflection components 
by measuring the ratio between the reflected and total ob- 
served flux. 

The blurred reflection model illustrated so far is able to re- 
produce all the spectral complexity of Ark 120 leaving no 
clear structure in the residuals (see Fig. |3). However, it is 
necessary to discuss in more detail the best-fitting param- 
eters that come to light (Table [2]). In particular, the disc 
emissivity profile is suggested to be very steep (q > 6.2, con- 
sistent with the largest value allowed in the model, which is 
10). This would imply that energy dissipation is extremely 
concentrated in the disc inner regions, virtually all the emis- 
sion arising within ~2-3 r g . Also, an inclination of 57° seems 
to be too large for a source with such a clean line of sight, 
although a misalignment of the accretion disc with the ob- 
scuring material is possible (see Lawrence & Elvis 2010)0 

3 We note that the host of this Seyfert nucleus is a low-inclination 
(i ~ 26°) spiral galaxy (Nordgren et al. 1995). 
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Figure 4. Contour plots (68, 90, 99 per cent confidence levels) 
showing the degeneracy among the critical blurring parameters. 
The emissivity index (horizontal axis) appears to be poorly con- 
strained against the disc inclination (vertical axis, upper panel) 
and inner radius (vertical axis, lower panel). The crosses mark 
the best-fitting solution when all the three parameters are free to 
vary. 

Indeed, it turns out that at this level of signal-to-noise the 
critical blurring parameters can be highly degenerate with 
each other, as shown in Fig. 3] Specifically, the emissivity 
index appears to be poorly constrained with respect to the 
inner radius and the disc inclination. We have therefore per- 
formed again all the fitting procedure, after freezing in turn 
each one of these three key parameters to a more convincing 
value, that is q — 5, n n = 3 r s and 9 = 40° (all obtained from 
the intermediate analytical steps). In all of these attempts 
we are able to recover a fit that is statistically equivalent to 
the one of reference: according to an F-test, the ancillary fits 
match the general case at a confidence level of 26, 15 and 31 
per cent, respectively^ In the light of these considerations, 
we can conclude that the blurred reflection model provides 
a valuable interpretation of the broadband X-ray emission 
of Ark 120 without requiring any extreme parameter. 
Owing to the degeneracy discussed above, it is not possible 
to put solid constraints on the black hole spin, as recently 
done e.g. in the case of Fairall 9 (Schmoll et al. 2009). In 
this effort we have switched to a different relativistic kernel 
(kerrconv; Brenneman & Reynolds 2006), which allows 
the dimensionless spin parameter a = cJ/GM 2 to vary, as- 
suming that the disc extends down to the innermost stable 
circular orbit (J and M are the black hole angular momen- 
tum and mass). The former q-r ln degeneracy is shifted into 
a q-a degeneracy; yet, by freezing again q = 5 we find that 
0.24 < a < 0.93 at a 99 per cent confidence level, with a 

4 The caveats against the use of the F-test for X-ray spectral 
analysis (e.g. Protassov et al. 2002) mainly concern the detection 
of marginal lines. This is not the case here, and the F-test is 
only intended to assess the degree of degeneracy among the key 
parameters. 



© 2010 RAS, MNRAS OOO.ITJITTI 



6 E. Nardini et al. 



best estimate of a ~ 0.74. It should be noted that the blur- 
ring parameters required to model properly the soft excess 
do not show an obvious agreement with those fitting sep- 
arately the broad iron line through the same LAOR profile, 
as pointed out in the previous section. In the framework 
of a pure disc reflection model, a substantial discrepancy in 
the blurring/ionization parameters of the soft excess and the 
companion broad iron K line is not expected. Anyway, we re- 
mark that the critical quantities listed in Table [2] should not 
be taken at their face value, due to the large degeneracy. In- 
deed, even after freezing ri n = 10 r g we still obtain a fully ac- 
ceptable fhflwith Xv — 0.89. In addition, it has already been 
suggested that the line detection can be difficult with present 
data quality unless some favourable conditions are met, e.g. 
a high iron abundance and/or a reflection-dominated spec- 
trum (as for 1H0707-495; Zoghbi et al. 2010). By inspecting 
carefully the residuals in the ~5-7.5 keV energy range, ten- 
tative evidence of some unfitted structure is found. We have 
therefore frozen the best-fitting model altogether and tried 
to add a gaussian profile, realizing that the residuals are sig- 
nificantly smoothed after the marginal detection of a narrow 
line centred at 6.69(±0.06) keV, whose EW is < 29 eV. If 
included in the model from the beginning, this line cannot 
be distinguished with such precision|f| Using Fig. [2] as a ref- 
erence, it seems plausible that the iron K-shell feature in 
Ark 120 consists of three narrow components (due to neu- 
tral, He-like and H-like iron) plus a faint relativistic disc line, 
whose parameters cannot actually be constrained against 
the blurring of the soft excess. 

It is significant that the narrow lines are not accounted 
for by the self-consistent reflection model we have adopted. 
This notwithstanding, an origin different than fluorescence 
in material exposed to intense X-ray illumination is diffi- 
cult to consider. We stress that only the putative H-like 
line keeps unchanged the properties (i.e. a central energy 
of ~6.97 keV and an EW of ~30 eV) already found in the 
analysis restricted to the iron region. Conversely, the line 
at lower energy is no longer fully consistent with neutral 
iron emission (it peaks at ~6.46 keV), and its EW is sig- 
nificantly reduced (~38 eV). This is because genuine neu- 
tral iron emission is now embodied in the cold reflection 
component (see Fig. [3]). These two lines then point to in- 
termediate/high ionization states of iron, and may arise in 
a much more internal environment. In order to probe this 
possible location, we have again left the width of these lines 
free, achieving a marginal improvement (Ax 2 — —3.7 with 
the loss of one degree of freedom, so that the probability 
of chance improvement is non- negligible). All the previous 
best-fitting parameters are remarkably confirmed, and the 
resulting width of a ~ 80(±50) eV corresponds to a FWHM 
velocity of 8.7(±4.3) x 10 3 km s _1 , which is consistent with 
an origin in (or slightly internal to) the optical BLR. The 
emission lines seem therefore to be connected to reflection 
neither at the pc nor at the r s scale. A possible explana- 
tion hints at some blobs of material orbiting in the outskirts 

5 Interestingly, in this case the disc inclination is ~25°, in excel- 
lent agreement with that of the host galaxy. The emissivity index 
is steep (~7.9) but again is poorly constrained. 

6 Moreover, even if a third narrow line is added to the best fit of 
the iron feature alone obtained before, the fit does not improve 
with any statistical significance. 



of the disc, forming a sort of X-ray BLR as suggested by 
Blustin & Fabian (2009), or where the height-radius rela- 
tion of the disc changes. 

Finally, we briefly comment on the relative amplitudes of the 
reflection components, which have been computed in terms 
of fractional contributions to the observed flux in both the 
XIS and the HXD energy range. Shortward of ~10 keV, the 
cold reflector accounts for ~4 per cent of the received emis- 
sion, its significance growing at higher energies where the 
Compton hump comprises ~18 per cent of the total flux 
(or equivalently ~30 per cent with respect to the power-law 
component only). This strength can be compared to that ex- 
pected from a slab subtending a solid angle of 2-7T at the illu- 
minating source by performing a straightforward test with a 
simple powerlaw+pexrav model (Magdziarz & Zdziarski 
1995) in xspec. Recalling that pexrav includes only the re- 
flected continuum, hence a sound comparison with reflionx 
is possible only beyond ~2 keV, we find a remarkable agree- 
ment between our measure and the 'dummy' estimate. It 
is worth noting that even in the presence of strong light 
bending effects in the inner regions of the disc, at large dis- 
tance from the X-ray source isotropy should be recovered: it 
is therefore desirable for the cold reflector identifiable with 
the obscuring torus to see the same illuminating radiation as 
the observer at infinity. Such a prescription is entirely com- 
plied with in this case. The fractional contribution of the 
warm reflection component from the disc is instead rather 
flat, increasing with energy from ~19 at 0.5-12 keV to only 
~22 per cent at 12-40 keV. Again, these values do not pose 
severe problems in terms of the required reflection efficiency 
below ~10 keV, due to the wealth of soft X-ray emission lines 
that are not implemented in the standard pexrav/pexriv 
models; besides this, it is reasonable to allow for moderate 
light bending involving the blurred component. 
In conclusion, a model consisting of both a warm reflec- 
tion component arising from the photoionized disc surface, 
blurred by relativistic effects, and a cold one ascribable to 
a much more distant reprocessor proves to be successful 
in reproducing the broadband X-ray spectrum of Ark 120. 
Limited to this peculiar source, this scenario stands as the 
most convincing among those presented so far, also due to 
the minimal set of geometrical and physical assumptions in- 
volved. 

3.3 Alternative models 

As we have already mentioned in the introduction, several 
other models have been adopted in the last decade to fully 
reproduce the high-quality ~0.5-10 keV spectra of AGN 
provided by Chandra and/or XMM-Newton, specifically con- 
cerning the soft excess component. Over the years, most of 
these interpretations have been either discarded as too phe- 
nomenological or disputed because of some physical limita- 
tions. Nevertheless, for the sake of a fitting comparison ex- 
tending up to ~40 keV against the blurred reflection model, 
we have also explored some of the most common alterna- 
tives, on which we give a brief report in this section. All the 
models listed below differ from that based on blurred reflec- 
tion only in the ingredients apt to describe the soft excess. 
In particular, they all require a cold reflection component 
with iron abundance ~1. 1—1.5 (±0.4), accounting for ~4-6 
and ~20-30 per cent of the XIS and HXD/PIN observed 
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Table 2. Best-fitting parameters for the blurred reflection model of the whole 0.5-40 keV Suzaku spectrum. 
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a Photon index. 

b Iron abundance, in solar units. 
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d Ionization parameter of the blurred reflector, in erg cm s — 1 . 
e Reflection fraction over the 0.5-10 keV range. 
* Observed 0.5-10 keV flux, in 10 — 11 erg cm~ 2 s" 1 . 



flux, respectively, and also for the narrow emission feature 
at 6.4 keV. Other two unresolved gaussian lines are always 
needed in the 6.47-6.70 and 6.93-7.01 keV range: even if 
loosely consistent with He-like iron emission, the former is 
apparently trying to compensate somehow for the lack of the 
broadened disc line. In the following, we just provide the key 
parameters of the soft excess component in some illustrative 
cases. 

1) Blackbody (xS/d.o.f. = 0.995/471): the best-fitting disc 
temperature is ~0.14 keV, and therefore lies in the range 
characteristic to the thermal scenario. No significant im- 
provement can be obtained by adopting a more realistic (e.g. 
multicolour) blackbody profile. 

2) Broken power law (xS/d.o.f. = 0.963/471): as one can 
immediately see from Fig. [T] a double-sloped power law is 
sufficient to achieve a good fit to the low-energy spectrum. 
The soft photon index is ~2.33, the hard one is ~2.00; the 
break-point energy is established at 1.73 keV. 

3) Smeared absorption (x^/d.o.f. = 1.02/470): we have used 
the SWlNDl model to mimic absorption in a partially ionized 
relativistic gas, even if the latest simulations of the velocity 
and density structure of any possible accretion disc wind or 
outflow rule out this scenario as the origin of a smooth soft 
excess (Schurch et al. 2009). The column denisty of the pu- 
tative absorber is Nn ~ 17 x 10 22 cm -2 , and the ionization 
parameter is log£ ~ 3.3. The gaussian velocity dispersion 
in this warm gas is the maximum allowed, i.e. v/c — 0.5. 
Since the power-law index is slightly steeper than usual, be- 
ing ~2.12, this model presents larger deviations from the 
data at the higher energies. 

4) Cold Comptonization (xl/d.ol. = 0.918/470): many 
codes with a different degree of complexity are available 
to describe Comptonization (e.g. Poutanen and Svensson 



1996). However, it is not possible to put firm constraints on 
the plasma properties and geometrical structure with the 
present data quality and energy coverage; we have there- 
fore adopted the essential COMPTT code of Titarchuk (1994). 
The exact temperature of seed photons is not very impor- 
tant, provided that its value is reasonably low, so it has 
been frozen to 50 eV as broadly presumable for a source 
like Ark 120. The cold corona turns out to have an electron 
temperature of ~0.28 keV with an optical depth of ~13, 
even if the two parameters are expected to be degenerate to 
some extent. This optically thick soft component has a lumi- 
nosity of ~10 43 erg s _1 , and can be roughly approximated 
by a blackbody. This gives a size of the emitting region of 
~4x 10 11 cm, whereas r g ~2x 10 13 cm (the mass of the 
black hole in Ark 120 is estimated from reverberation map- 
ping to be M BH ^ 1.5 x W s Mq; Peterson et al. 2004). 
Summarizing, the Comptonization scenario can be regarded 
at present as the only physical alternative to the blurred 
reflection model in order to account for the soft excess of 
Ark 120, but its details are difficult to probe; the argued 
compactness of the soft X-ray emitting region represents 
a further problem for the nature and the geometry of the 
Comptonizing plasma. 

3.4 Timing analysis 

In addition to spectral fitting, the other fundamental ap- 
proach to understand the properties of the central engine 
in AGN is the study of variability. In principle, the latter 
can also provide independent information to discriminate 
between the cold Comptonization and the blurred reflection 
scenario for the soft excess, since the two models make quite 
divergent predictions about the behaviour of the source in 
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the time domain. In a very simplistic way, the time lags af- 
fecting the different energy bands that have been observed 
in many AGN can be regarded as a natural consequence of 
a Comptonization process: the time delay between the hard 
and the soft X-ray variations is due to the larger number 
of scattering events that the high-energy photons have to 
experience before escaping to infinity. However, matters ap- 
pear to be much more complicated than this naive picture 
(e.g. Arevalo & Uttley 2006; McHardy et al. 2007). On the 
other hand, if the soft excess is a reflection signature, the 
soft band should lag behind the hard power-law component 
as expected in a reverberation context, after singling out the 
timescale corresponding to the light crossing time from the 
primary source to the reflector and back. Evidence in this 
sense has been recently found in 1H0707-495 (Fabian et al. 
2009; Zoghbi et al. 2010). 

Although Suzaku is not an ideal observatory to carry out 
detailed timing studies of AGN, mainly because of the fre- 
quent gaps induced by the short orbital period, it is worth 
checking the variability pattern of Ark 120 for possible hints. 
The total XIS light curve is shown in Fig. and displays a 
very gentle variation, in part expected due to the large black 
hole mass. Using a time resolution of 500 s, the fractional 
rms variability amplitude F vilI (defined as in Vaughan et al. 
2003) is 9.1 (±0.2) per cent over the 0.5-10 keV ranged In 
order to check for a possible energy dependence, we have 
computed the individual variability amplitudes in four en- 
ergy bands (0.5-1, 1-2, 2-5 and 5-10 keV), which turn out 
to be 8.7(±0.4), 9.7(±0.2), 9.1(±0.3) and 8.5(±0.7) per cent, 
respectively. This does not provide any straightforward indi- 
cation for different variability patterns, and these are likely 
to be actually consistent with each other. Moreover, we note 
that the four bands above cannot be clearly linked to dif- 
ferent spectral components, since the direct power law is 
dominant over the whole spectrum of Ark 120. Incidentally, 
the shape of the HXD/PIN light curve closely resembles 
the trend observed in the XIS range, although in this case 
Jvar = 19.1(±3.4), due to the larger noise. A similar indica- 
tion in favour of a coarse coherence with no lags among the 
energy bands comes from the 1-2, 2-5 and 5-10 over 0.5- 
1 keV hardness ratios (Fig. [5]), which are all well- fitted with 
a constant value (x 2 /d.o.f ~ 246/231, 252/231 and 267/231, 
corresponding to null hypothesis probabilities of 24, 16 and 
5 per cent, respectively). The lack of differential variability 
was also found by De Marco et al. (2009), who have analyzed 
the XMM-Newton spectrum focusing on the iron K line re- 
gion. The total time elapsed in this Suzaku observation is 
~185 ks, corresponding to ~1.2 keplerian orbital periods at 
10 r g (e.g. Bardeen, Press & Teukolsky 1972). It is therefore 
unlikely to reveal dramatic changes on such timescales, de- 
spite the small radii involved. 

In short, the timing analysis of Ark 120 is inconclusive and 
does not point to any specific physical model. Yet, by con- 
firming that the source variations involve the flux amplitude 
only, it offers substantial support to our previous spectral 



7 All the values reported in this section have little dependence 
on the adopted time bin. This does not undermine the general 
conclusions: in this source, the light crossing time over a distance 
of 1 r g is ~750 s, so that most of the variability is expected to 
take place on larger timescales. 
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Figure 5. Top panel: 0.5-10 keV background-subtracted light 
curve of Ark 120, summed over the three XIS detectors and plot- 
ted adopting a time resolution of 500 s. The flux variation is 
rather gentle, with an overall rms variability amplitude of ~9 per 
cent. Lower panels: respectively from above, 1-2, 2-5 and 5—10 
over 0.5-1 keV hardness ratios, which turn out to be remarkably 
stable over the entire observation; the red solid lines show the 
best-fitting constant values. (Note: the outlying points with large 
error bars correspond to the ends of a single Suzaku orbit, and 
are just artifacts of the time binning). 

study, given that the average X-ray emission genuinely rep- 
resents the single physical state in which the source was 
caught during the Suzaku observation. 

3.5 The 2003 XMM-Newton observation 

Previously to Suzaku, Ark 120 was monitored by XMM- 
Newton on 2003 August 24-25. A detailed study of this 
observation is presented in Vaughan et al. (2004). In particu- 
lar, a wide range of models (including reflection components, 
disc blackbodies, bremsstrahlung and thermal Comptoniza- 
tion) were tested in order to give an adequate description of 
the broadband spectrum; yet, none of these was able to ac- 
count for the properties of the soft excess. It is not our aim to 
perform the same comprehensive analysis over again, repro- 
ducing equivalent results. Since then, anyway, the available 
reflection models have been substantially upgraded, after 
Ross & Fabian (2005). So it is worth reviewing the very high- 
quality XMM-Newton observation in the light of the blurred 
reflection scenario we have developed in this work. Our re- 
duction follows that illustrated in Vaughan et al. (2004), and 
yields to a full agreement of the final data products. The fol- 
lowing analysis concerns the EPIC-pn spectrum only, and is 
confined to the 0.5-10 keV energy range for consistency with 
the Suzaku XIS operational window. Moreover, due to the 
large count rate in the soft band the ionization and blurring 
parameters are highly sensitive to calibration uncertainties, 
which can be quite large below ~0.6 keV (see also Papadakis 
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Table 3. Best-fitting reflection parameters for the 2003 XMM-Newton observation. 



r 


4-0 01 1 
Z.11D_ oog 


£b,l 


ion 
299Iig 






?b,2 




E 2 


6 64+ 03 
u - u ^-0.04 




r 9+0.4 


E 3 


7.00 ± 0.03 


Hn 


2 26+ - 11 


17 


10* 




50* 


EWi 


47+23 
^'-22 




< 0.008 


EW 2 


26 ± 7 


Rb,l 


n 91H +0.013 


EW 3 


25 ± 15 


Rb,2 


0-062™ 


Fe 


0.75* 




6.80t°;« 


€d 


1.0* 


xl 


2044/1857 



et al. 2010). 

The spectral shape and flux level in the two observations 
are remarkably similar, hence we have directly applied the 
best-fitting model obtained for the Suzaku case. Iron abun- 
dance and disc inclination are obviously not time-dependent 
parameters, and have been frozen to 0.75 (solar units) and 
50° , respectively. Under these assumptions we obtain a rea- 
sonable fit to the data, even if not formally acceptable in 
a statistical sense (xt ~ 1.1G), as quite usual at the high- 
est levels of data quality. A further improvement can be 
achieved by allowing for a slightly larger degree of complex- 
ity in the model. First, a third narrow line is included for 
He-like iron emission as also pointed out by Brenneman & 
Reynolds (2009). Second, we introduce an additional reflec- 
tion component, subject to the same relativistic smearing 
but with an independent ionization parameter, in order to 
allow for possible non-uniformity within the disc. As a re- 
sult, this extra component turns out to require a low ion- 
ization, while the cold distant reflection becomes negligible. 
The value of xi/d.o.i. drops to ~1. 10/1857 and the best- 
fitting key parameters are very well-matched to those ob- 
tained for the Suzaku observation (see Table [3}. Conversely, 
beyond the iron feature there remains some excess emission 
that cannot be fitted by increasing the model complexity. 
These deviations of the model from the data are of the or- 
der of ~10 per cent at most (Fig. [6]). We note that if the disc 
inclination is left free to vary, a value around ~80° (see also 
Nandra et al. 2007) is able to offset any ~8-10 keV excess, 
while the iron abundance is not a critical parameter in this 
sense. Interestingly, the Galactic column density converges 
towards the upper limit of the admitted range: if instead Ah 
is frozen to the lower value the fit quality marginally wors- 
ens, but the free parameters smoothly resettle so that the 
picture above is not modified. We have also checked the be- 
haviour of the model down to 0.3 keV: a very clean bump in 
the data/model ratio emerges at 0.3-0.5 keV. After freezing 
the best-fitting model, this feature is quite well-reproduced 
by a blackbody of temperature ~40 eV, even if the residuals 
in the 0.3-0.7 keV range display some tentative structure. 
Finally, we have tested for completeness the cold Comp- 
tonization model: in this case, the parameters of the soft 
corona are in excellent agreement with those found above for 
the Suzaku spectrum, but the fit is less accurate (xl ~ 1.19) 
and large residuals remain in the 0.5-0.7 and 5.5-7 keV 
range. Our revision of the high-quality XMM-Newton ob- 
servation seems therefore to confirm the reliability of the 
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Figure 6. Top panel: as Fig. 0] but for the 2003 XMM-Newton 
observation, with the same code. The additional cold/blurred re- 
flection component is plotted in orange (short-dashed line). The 
inner box shows the profile of the iron feature, obtained as the 
ratio over a simple power-law trend. Bottom panel: data/model 
ratio to the best fit, emphasizing the slight excess emission long- 
ward of the iron feature. Errors bars are omitted for clarity. 



blurred reflection interpretation for soft X-ray emission of 
Ark 120. 



4 CONCLUSIONS 

We have reported on a long Suzaku observation of Ark 120, 
one of the few known examples of nearly absorption- free 
active nuclei. This 'bare' BLS1 galaxy is the optimal can- 
didate to investigate the nature of the ubiquitous soft ex- 
cess observed in the X-ray spectra of unobscured AGN. Its 
broadband 0.5-40 keV spectrum reveals a significant cur- 
vature with respect to the main power-law component, as 
excess emission is found at both low (< 2 keV) and high 
(> 10 keV) energies. Also, a complex iron K-shell feature is 
clearly detected at ~5-7 keV, showing a structured profile 
in which can be identified two (three, tentatively) narrow 
gaussian cores originating roughly at the BLR scale, and a 
relativistically skewed line arising from the disc at r > 10 r g . 
The latter, even if not as prominent as in the most impressive 
cases, cannot be interpreted through a blend of unresolved 
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lines or a system of multiple absorbers. Hence relativistic 
effects have to be included in the broadband analysis. 
We have shown that all the spectral complexity can be suc- 
cessfully explained in terms of a self-consistent reflection 
model, allowing for both a warm/blurred and a cold/distant 
reflection component. In this picture, the reprocessing of the 
primary X-ray radiation takes place in two distinct regions, 
that can be plainly identified with the inner accretion disc 
(r < 10 r g ) and the far-off obscuring torus. Depending on 
the exact geometry, however, additional physical locations 
(e.g. a belt of orbiting clouds) may contribute, since fluores- 
cent iron emission suggests the existence of a wide range of 
ionization states. This blurred reflection scenario is further 
corroborated after reviewing a previous high-quality XMM- 
Newton observation of Ark 120 . Therefore, as for this single 
source, the present interpretation stands as the most con- 
vincing among those that have been proposed so far, due 
to the minimal set of geometrical and physical assumptions 
involved. 

In general it is difficult to rely on comparably clean lines 
of sight to the nuclear regions, and we expect that other 
processes (likely non-thermal Comptonization, warm ab- 
sorption, reprocessing in winds/outflows) also play a non- 
negligible role in shaping the X-ray emission of type 1 AGN. 
As a consequence, our aim in the near future is to extend this 
study to a large sample of low-obscuration sources, possibly 
observed by both Suzaku and XMM-Newton to combine the 
high-energy coverage of the former and the large effective 
area of the latter. This will allow us to infer the average 
contribution of blurred disc reflection to the soft excess and 
to the X-ray luminosity of AGN, and at the same time to 
test the basic predictions of the light bending model and to 
explore the nature of the illuminating source itself. We note 
that if blurred reflection is really the dominant component 
of the soft excess, the black hole spin distribution may be 
biased towards the larger values; indeed, it has already been 
suggested in many previous works (e.g. Volonteri et al. 2005) 
that the supermassive black holes in AGN should be rapidly 
rotating. Even if the spin cannot be firmly constrained in 
Ark 120 due to the large degeneracy among the blurring 
parameters (but provisionally a ~ 0.7), our case study is 
apparently consistent with this scenario. 
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